Probing low- a? QCD with cosmic neutrinos at the Pierre Auger Observatory 
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The sources of the observed ultra-high energy cosmic rays must also generate ultra-high energy 
neutrinos. Deep inelastic scattering of these neutrinos with nucleons on Earth probe center-of-mass 
energies yfs ~ 100 TeV, well beyond those attainable at terrestrial colliders. By comparing the rates 
for two classes of observable events, any departure from the benchmark (unscreened perturbative 
QCD) neutrino-nucleon cross-section can be constrained. Using the projected sensitivity of the 
Pierre Auger Observatory to quasi-horizontal showers and Earth-skimming tau neutrinos, we show 
that a 'Super-Auger' detector can thus provide an unique probe of strong interaction dynamics. 



A new "window on the Universe is expected to open 
soon through the observation of ultra-high energy cosmic 
neutrinos by detectors such as the Pierre Auger Observa- 
tory and IceCube Q ■ The study of their interactions 
in the Earth's atmosphere and crust is particularly inter- 
esting as this provides a probe of QCD in the kinematic 
region of very small values of Bjorken- x: x < 10~ 3 , where 
conventional calculations done in the DGLAP framework 
become inadequate. 

The surprising discovery at HERA that the gluon dis- 
tribution function rises sharply with decreasing x for 
high virtuality of the exchanged gauge boson Q 2 , im- 
plies a strong increase in the neutrino-nucleon total cross- 
section. This would however violate unitarity (the Frois- 
sart bound) if continued indefinitely, hence there must be 
a departure from this behavior at very low values of x. 

Various possibilities have been entertained in this con- 
text (see e.g., :3j for a recent review) which predict a 
slower rise of the cross-section than the usual DGLAP 
evolution based calculation £| commonly adopted for es- 
timating the effective aperture of detectors, e.g. of Auger, 
which is sensitive to neutrinos of energy > 10 s GeV |5|. 
Thus, a measurement of the neutrino-nucleon cross- 
section at such energies, even if relatively crude, would 
nonetheless be of great interest. For example one can 
directly test models where saturation effects cause the 
gluon distribution function to freeze or even decrease with 
x below some threshold, x < 1CU 5 . The former implies a 
v — N cross-section which stays about a factor of 2 be- 
low the standard expectation for E v ~ 10 8 — 10 10 GeV, 
where Auger is most sensitive. In the latter case the 
cross-section remains sensibly constant with increasing 
energy above ~ 10 9 GeV, so has a value 10 times below 
the usual expectation at the upper end of this range. 

One might wonder how such a measurement can be 
done given the large uncertainties in the expected fluxes 
of (the yet to be detected!) cosmic neutrinos. The op- 
portunity arises because of the discovery that the muon 
and tau neutrinos are maximally mixed. Hence cos- 
mic neutrino beams, generally expected to be z^'s and 
i/ e 's from the decays of pions, kaons and perhaps heavy 



flavours, will necessarily contain a substantial compo- 
nent of 2/ T 's by the time they reach Earth [||. It has 
been noted that Earth-skimming i/ T 's will gener- 

ate upward going air showers when they interact in the 
crust. By contrast neutrinos of all flavours will gener- 
ate deeply penetrating quasi-horizontal air showers which 
are distinctive in having an electromagnetic component 
unlike hadron-initiated air showers at such large incli- 
nations 01 • The crucial observation is that while 
the rate of quasi-horizontal z^-showers is proportional to 
the v — N cross-section, the rate of detectable Earth- 
skimming z^-'s is not. This is because of several effects 
which come into play in the latter case. For example, 
a decrease in the cross-section will increase the number 
of r's produced in a region close enough to Earth's sur- 
face that are likely to escape. However, the resulting r's 
emerge preferentially at angles outside the acceptance 
typical of surface arrays. Consequently, the number of 
detected r's is relatively independent of the cross-section. 
All in all, the ratio of the two classes of events provides a 
measure of the absolute v — N cross-section, even when 
there are large uncertainties in the incoming cosmic flux. 

For example although the (so far unknown) sources of 
the observed ultra-high energy cosmic rays must also be 
sources of ultra-high energy neutrinos, it is difficult to 
calculate the expected neutrino flux in terms of the ob- 
served cosmic ray flux, given our ignorance of the opacity 
of the sources. The usual benchmark here is the so-called 
(all flavours) Waxman-Bahcall flux 
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for 10 3 < E v /QeV < 10 11 , derived assuming that the 
sources are 'transparent' and that 60% of the proton en- 
ergy is converted to pions. We will use this flux to esti- 
mate the event rates for various models of the neutrino 
scattering cross-section. 

The cross-section for charged current (CC) v — N scat- 
tering is |l3| 
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is the differential cross-section given in terms of the struc- 
ture functions F 2 , Fl and a;F 3 , and V + = 1 + (1 — y) 2 , 
Y- = 1 — (1 - y) 2 with y = Q 2 /sx (s = 2E u m N is the 
center-of-mass (cm.) energy, Gf is the Fermi constant 
and mw is the VF-boson mass). For simplicity, we are 
considering only CC interactions, as neutral current (NC) 
interactions are subdominant. At leading order (LO) in 
perturbative QCD, the structure functions are given in 
terms of parton distributions as F 2 = J2i x [q%( x i Q 2 ) + 
qi(x,Q 2 )}, xF 3 = Y^i x [Qi{ x ,Q 2 ) + q~i( x ,Q 2 )} and F L = 
0. However at NLO these relationships involve further 
QCD-calculable coefficient functions and contributions 
from F^ can no longer be neglected. Parton distribu- 
tion functions (PDFs) are determined in fits to deep in- 
elastic scattering (DIS) data by the following procedure. 
The parton distribution functions are parameterised at 
some initial scale Qo ~ 1 GeV and then evolved, using 
the next-to-leading order (NLO) DGLAP equations, to 
higher values of Q 2 ; they are then convoluted with QCD- 
calculable coefficient functions to give NLO predictions 
for the structure functions, which are then fitted to the 
DIS data Such fits have been made by several dif- 
ferent groups 0, 0, 0| ; recent analyses have included 
estimates of the uncertainties on the PDFs coming from 
experimental uncertainties. 

As Q 2 increases, the parton distribution functions 
(particularly of the gluon) grow due to QCD evolution, 
so that the neutrino cross-section will also grow until the 
propagator cuts off the growth at Q 2 ~ m |/. Hence the 
typical x value probed is x ~ j2m,^E v |l7|. For neu- 
trino energies E v ~ 10 8 — 10 10 GeV this translates into 
small x values of 1(T 4 - 1(T 6 at Q 2 ~ 10 4 GeV 2 . HERA 
measurements do extend down to i ~ 10~ 6 , but only 
at Q 2 < 0.1 GeV 2 , while for Q 2 - m 2 w the LHC will 
probe 5 x 1CP 4 <i<5x 10~ 2 . To probe down to the 
same kinematic region as Auger would require a hadron 
collider with cm. energy exceeding 10 3 TeV. 

At small x and high Q 2 the v — N cross-section is 
dominated by sea quarks produced by gluon splitting 
g — ► qq. In this kinematic region, the parametrisation 
of the gluon momentum distribution is approximately: 
xg(x, Ql) oc x~ x for x < 1, where A ~ 0.3 - 0.4. The re- 
sulting CC v — N cross-section was originally calculated 
at leading order using 1996 parton distribution functions 
and parametrised as: tr„° r = 5.53 (£^/GeV) a363 pb for 
10 7 < E v < 10 12 GeV Hf. This is the benchmark un- 
screened cross-section widely used to evaluate sensitivi- 
ties of ultra-high energy cosmic neutrino detectors. We 
show this in Fig. ^ along with our updated calculation 
using a modern PDF analysis |16| which included the 
final data from the first phase of HERA running (1993- 
2000); this data is essential to obtain information on the 



PDFs at low- a;. Our NLO analysis includes corrections 
for heavy flavours, and, most importantly, a full treat- 
ment of experimental uncertainties which were not con- 
sidered in [j]. We also show in Fig.^a simple power-law 
fit to our improved CC v — N cross-section: 

a* L ° r = 6.04 ± 0.40 (^/GeV)°- 358±0 005 pb . 



(4) 



for 10 7 < E v < 10 12 GeV. (The tabulated cross-section 
and other details are provided elsewhere [T^. - ) 

However, when x is small two further considerations 
are relevant. First, when x is sufficiently small that 
a s ln(l/a;) ~ 1, it is also necessary to resum these 
large logarithms, using the BFKL formalism. Second, 
the gluon density at small- a; is very high and both the 
DGLAP and the BFKL formalisms neglect non-linear 
screening effects due to gluon recombination. Such effects 
would tame the rise of the gluon distribution function at 
small x and may even lead to saturation. An efficient 
way of modelling this is the color dipole framework in 
which DIS at low x is viewed as the interaction of the qq 
dipole to which the gauge bosons fluctuate. An unified 
BFKL/DGLAP calculation supplemented by screening 
effects, as well as nuclear shadowing (following the cal- 
culation in [19j for A = 12), predicts a decrease of the 
cross-section by a factor of ~ 1.2 — 2 in the relevant 
energy range [2(|. An alternative recent approach uses 
the colour glass condensate formalism [2l]]; this predicts 
a similar suppression when a dipole model [22| which 
fits data from RHIC is used |2^|. The predicted cross- 
section tr™ is even lower if a different dipole model de- 
veloped to fit the HERA data [24[ is used and the gluon 
distribution is assumed (more speculatively) to decrease 
for x < 10~ 5 |2^| . As seen in Fig. ^thrs is a factor of 
^2 — 20 below the unscreened cross-section in the rele- 
vant energy range. Other possibilities for the high energy 
v — N cross-section have been discussed [25j ; an exciting 
development formulates DIS using gauge/string duality 
and provides new insights into low- a; dynamics |26|. 

The Pierre Auger Observatory is the largest cosmic 
ray detector in the world [J and when complete will oc- 
cupy two sites — one in the Southern and one in the 
Northern hemisphere. Construction of Auger South on 
a plateau in Western Argentina is well advanced and it 
is operating in hybrid mode employing fluorescence de- 
tectors overlooking a ~ 3000 km 2 ground array of water 
Cerenkov detectors. In addition to studying the highest 
energy cosmic rays, Auger is also capable of observing 
ultra-high energy cosmic neutrinos 0, Q • 

Given an isotropic flux of neutrinos W the rate of 
quasi-horizontal (QH) showers expected to be observed 
at Auger is proportional to the v — N cross-section: 



A/qh oc J dE sh g vN {E u ) A QH {E sh ,9 z ) <j> v (E v ) 



(5) 



Here Aq^E^, 6 z ) is the Auger acceptance, which de- 
pends on the zenith angle 9 Z |27|. and E^ = yE v for 
and j/ t , whereas E^ = E v for v e . The expected event 
rate for the WB flux (EqQ is given in Table [I] 
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Neutrino Energy (GeV) Number of Earth— skimming events 



FIG. 1: Predicted cross-sections for neutrino- nucleon scatter- 
ing at high energies. The line with its la error band is the fit 
(Eq0J to our calculated cr^ s °(-E„) (points with error bars). 
For comparison we show a^2s(Ev) (dotted line) Bl, a ^(E v ) 
(dashed line) H|, and ofc X {E v ) (dot-dashed Une"|23||. 



TABLE I: Neutrino events per year at Auger for different 
models of the v — N cross-section, adopting two benchmark 
cosmic fluxes and shower energy threshold of I0 8 (f0 9 ) GeV. 



Model 


Waxman-Bahcall 


Cosmogenic 




ATqh 


A/es 


A/" QH 


A/es 


LO 
^unscr 


0.15 (0.092) 


3.0 (0.62) 


0.061 (0.039) 


1.2 (0.35) 


NLO 
^unscr 


0.14 (0.080) 


3.0 (0.61) 


0.057 (0.036) 


1.2 (0.34) 


KK 

^scr 


0.10 (0.057) 


2.7 (0.54) 


0.042 (0.027) 


1.1 (0.31) 


HJ 

® scr 


0.048 (0.022) 


1.8 (0.32) 


0.018 (0.010) 


0.7 (0.18) 



The situation is different for showers initiated by r's 
created by CC interactions of Earth-skimming (ES) ^ r 's. 
To a first approximation the number of such events is 

A/es oc J dE sh dcosf? d$ P{9,cj ) )A E s{E sh ,9)^{E u ) , 

(6) 

where 

p(M)= r*F e ~* /r e [*-^-M] 

is the probability for a v T with incident nadir angle 9 
and azimuthal angle <f> to emerge as a detectable r. Here 
l T ~ 10 km is the typical r path length Q, I = 2i?e cos 9 
is the chord length of the intersection of the neutrino tra- 
jectory with the Earth (of radius ss 6371 km), 
and /£, ot are the CC and total neutrino mean free paths, 
respectively, and AEs(-/?shj$) is the experimental accep- 
tance which has a strong dependence on the angle, since 



FIG. 2: The expected number of Earth-skimming and quasi- 
horizontal neutrino events above 10 s GeV for different mod- 
els of the v — N cross-section; for each model, one of the 
lines assumes the cosmogenic spectrum and the other line 
the Waxman-Bahcall spectrum. We show as squares and 
stars respectively, the corresponding hypothetical measure- 
ments (with la statistical errors) that could be made in 10 yr 
with an Auger-like detector scaled up to 10000 square miles. 



the surface detector array can only see events within a 
few degrees of the horizon [lfj. In fact, the analytic ex- 
pression above is an oversimplification; it does not allow 
for r regeneration in the Earth and the r path length 
is not really a step-function. To take such details into 
account, we have carried out a simple Monte Carlo sim- 
ulation assuming the NC cross-section to be 40% of the 
CC cross-section for all models. The resulting ES event 
rates are given in Tabled- we find good agreement with 
results from a sophisticated Monte Carlo that models the 
environment of Auger and its acceptance accurately [28| . 

To evaluate the sensitivity to the assumed spectral in- 
dex we also consider the "guaranteed" cosmogenic neu- 
trino flux which has a peaked distribution in the energy 
range of interest [2^| . As shown in Table HJ if we con- 
sider events with E^ > 10 8 GeV, the change in the 
spectrum produces a variation in the ratio A/qh/A/es 
of less than 5%. We have also verified that a more 
steeply falling flux oc E~ 2 - 54 (with the same BO) causes 
a change in A/qh/A/es by about 10% (e.g, for cr^°, 
A/qh = 0.69 yr- 1 and A/es = 12 yr" 1 ). Such a flux 
is expected |3JJ if extragalactic cosmic rays from 'trans- 
parent' sources begin dominating the observed spectrum 
at 10 9 - 6 GeV rather than at ~ 10 10 5 GeV as is usu- 
ally assumed (see Fig. 5 in [22j for a comparison of these 
fluxes). Thus we conclude that the ratio A/"oh/A/"f,r pro- 
vides a robust estimate of the v — N cross-section 321. 
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Now we discuss the precision with which such a mea- 
surement can be made. In Fig. [21 we show the numbers 
of QH and ES events expected for different cross-section 
predictions and flux expectations. The data points are 
hypothetical measurements made over 10 years with a 
proposed 'Super Auger' array of area 10000 square miles 
(25600 km 2 ). Note that by observing just a dozen QH 
events, we can begin to distinguish between the theoreti- 
cal models, however from Table[I]we see that the number 
of events for the WB flux is an order of magnitude smaller 
in Auger even after 10 years. Hence we are led to enter- 
tain the idea of scaling up the array by an order of mag- 
nitude, perhaps by using radio detection methods [33|. 
Another possibility is to use satellite-borne fluorescence 
detectors such as EUSO and OWL which may attain the 
required order of magnitude increase in sensitivity |34| . 

The cosmic neutrino flux can be much higher than the 
conservative benchmark values we have adopted above. 
Data from HiRes |3lJ suggest that it may be [3(| just 
below the current experimental bound from AMANDA- 
BIO (35J. The expected event rates in Auger itself would 
then be high enough to test for any suppression of the 



UHE v — N cross-section. Since the spectrum in this 
case is softer than the WB flux, IceCube Q, with its 
lower energy threshold, should test this model very soon. 
Indeed IceCube has the sensitivity to see the benchmark 
WB flux within a few years. Thus there is an emerg- 
ing synergy between large cosmic ray arrays and cosmic 
neutrino detectors which will soon establish if there is a 
realistic possibility for exploring fundamental physics us- 
ing Nature's own high energy beams. Given that there 
is no conceivable terrestrial accelerator which can attain 
such energies, we believe that the construction of Super 
Auger should be considered seriously; the detection tech- 
niques are well developed and the resources required are 
no larger than for a contemporary collider experiment. 
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